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Abstract 



It is known that a boundary slip velocity starts to play important role when the length scale over 
which the fluid velocity changes approaches the slip length, i.e. when the fluid is highly confined, 
for example, fluid flow through porous rock or blood vessel capillaries. Zhu & Granick [Phys. Rev. 
Lett. 87, 096105 (2001)] have recently experimentally established existence of a boundary slip in a 
Newtonian liquid. They reported typical values of the slip length of the order of few micro-meters. 
In this light, the effect of introduction of the boundary slip into the theory of propagation of elastic 
waves in a fluid-saturated porous medium formulated by Biot [J. Acoust. Soc. Am., 28, 179-191 
(1956)] is investigated. Namely, the effect of introduction of boundary slip upon the function F(n) 
that measures the deviation from Poiseuille flow friction as a function of frequency parameter k is 
studied. By postulating phenomenological dependence of the slip velocity upon frequency, notable 
deviations in the domain of intermediate frequencies, in the behavior of F(k) are introduced with 
the incorporation of the boundary slip into the model. It is known that F(k) crucially enters 
Biot's equations, which describe dynamics of fluid-saturated porous solid. Thus, consequences of 
the non-zero boundary slip by calculating the phase velocities and attenuation coefficients of the 
both rotational and dilatational waves with the variation of frequency are investigated. The new 
model should allow to fit the experimental seismic data in circumstances when Biot's theory fails, 
as the introduction of phenomenological dependence of the slip velocity upon frequency, which is 
based on robust physical arguments, adds an additional degree of freedom to the model. If fact, it 
predicts higher than the Biot's theory values of attenuation coefficients of the both rotational and 
dilatational waves in the intermediate frequency domain, which is in qualitative agreement with 
the experimental data. Therefore, the introduction of the boundary slip yields three-fold benefits: 
(A) Better agreement of theory with an experimental data since the parametric space of the model 
is larger (includes effects of boundary slip); (B) Possibility to identify types of porous medium and 
physical situations where boundary slip is important; (C) Constrain model parameters that are 
related to the boundary slip. 

PACS numbers: 43.20.Jr: 47.27.Lx; 47.55.Mh; 47.60.+i; 62.65.+k 
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I. INTRODUCTION 



It has been a common practice in the fluid dynamics to assume that when a fluid flows 
over an interface with a solid, the fluid molecules adjacent to the solid have zero relative 
velocity with the respect to solid. So far, this widely used assumption, known as "no-slip 
boundary condition", has been successfully applied to the theoretical modeling of almost 
all macroscopic experiments. As relevantly noticed by Craig et al. the success of this 
assumption does not reflect its accuracy, but rather insensitivity of the experiment to a 
partial-slip boundary condition. It is known that the boundary slip becomes important only 
when the length scale over which the fluid velocity changes approaches the slip length, that 
is the distance behind the interface at which the fluid velocity extrapolates to zero, i.e. when 
the fluid becomes highly confined, e.g., blood flow through capillaries or fluid flow through 
natural porous rock. Recently, authors of Refs. |], || presented a convincing experimental 
evidence of a boundary slip in a Newtonian liquid. They performed direct measurements of 
the hydrodynamic drainage forces, which show a clear evidence of boundary slip. Also, they 
found that the boundary slip is a function of the fluid viscosity and the shear rate. These 
results have important implications for the blood dynamics in capillaries, the permeability 
of porous media, and lubrication of nano-machines. For example, results of Craig et al. 
suggest that red blood cells squeeze through narrow capillary walls more easily and induce 
less shear stress on capillary walls due to the boundary slip. Also, in oil production industry, 
the residual oil is difficult to produce due to its naturally low mobility. Thus, the enhanced 
oil recovery operations are used to increase production. It has been experimentally proven 
and theoretically validated that there is a substantial increase in the net fluid flow through 
a porous medium if the latter is treated with elastic waves 0, §], |5], f|. We may conjecture 
that the elastic waves via the pore wall vibration cause boundary slip of the residual oil 
droplets, which likewise red blood cells, squeeze through pores with less resistance, effectively 
increasing permeability of the porous medium. 

A quantitative theory of propagation of elastic waves in a fluid-saturated porous solid 
has been formulated in the classic paper by Biot M. After its appearance, this theory has 
seen numerous modifications and generalizations. One of the major findings of Biot's work 
was that there was a breakdown in Poiseuille flow above a certain characteristic frequency 
specific to this fluid-saturated porous material. Biot theoretically studied this breakdown 
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by considering the flow of a viscous fluid in a tube with longitudinally oscillating walls under 
an oscillatory pressure gradient. Biot's theory can be used to describe interaction of fluid- 
saturated porous solid with sound for a classic Newtonian fluid assuming no-slip boundary 
condition at the pore walls holds. However, in the light of recent experimental results of Ref. 
0, revision of the classic theory is needed in order to investigate novelties bought about by 
the boundary slip. 

Biot's theory has been a successful tool for interpreting the experimental data for decades, 
but there are circumstances when it fails. Gist || performed ultrasonic velocity and attenu- 
ation measurements in sandstones with a variety of saturating fluids and compared the data 
with the predictions of the Biot's theory. Velocity data show systematic deviations from Biot 
theory as a function of pore fluid viscosity. Ultrasonic attenuation is much larger than the 
Biot prediction and for several sandstones is nearly constant for a three decade variation in 
viscosity. This behavior contrasts with synthetic porous media such as sintered glass beads, 
where Biot theory provides accurate predictions of both velocity and attenuation. The simi- 
lar claim has been also made by earlier experimental studies e.g. Ref. M. Even Biot himself 



acknowledged [10] that the model of a purely elastic solid matrix saturated with a viscous 
fluid is rather idealistic and there are rare occasions where it gives satisfactory agreement 
with the experiment. Thus, he generalized theory applying the correspondence principle, and 
incorporated the solid and fluid attenuation, and dissipation due to non-connected pores. 
In this work Biot introduced the viscodynamic operator, which is a sum of the viscosity 
and inertial terms, and it accounts for the dynamic properties of the fluid motions in pores 
valid for the both low and high frequencies. The theory can model other, more complex, 
mechanisms such as thermoelastic dissipation, which produce a relaxation spectrum, or local 
flow mechanisms as the 'squirt' flow fIT||, [fI2|. A time domain formulation based on mem- 



ory variables has been also successfully introduced |13[ , |T4| . In order to narrow down the 
source of discrepancy between Biot theory and experiment, Gist || tried to modify natural 
sandstone by curing a residual saturation of epoxy in the pore space, filling small pores 
and micro-cracks. This altered rock has a significantly reduced attenuation, demonstrating 
the dominance of small pores in controlling ultrasonic attenuation. Then he suggested that 
two non-Biot viscous attenuation mechanisms are needed: local flow in micro-cracks along 
grain contacts, and attenuation from pore-wall surface roughness. The model model based 
on these two mechanisms seemed to succeed in improvement of fitting the data. However, 
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this by no means does not exclude other possibilities, and all plausible effects that could 
explain the discrepancy should be explored. At low frequencies, the fluid flow is of Poiseuille 
type and the inertia effects are obviously negligible in comparison with the effects of the 
viscosity. At high frequencies, however, these effects are confined to a thin boundary layer in 
the vicinity of the pore walls and the inertia forces are dominant. Therefore, we conjecture 
that in such situations, apart from other yet unknown effects, non-zero boundary slip effect 
maybe responsible for the deviations between the theory and experiment. This justifies our 
aim to formulate a model that would account for the boundary slip. 

In the section II we formulate theoretical basis of our model and in section III we present 
our numerical results. In the section IV we conclude with the discussion of the results. 



II. THE MODEL 



In our model we study a Newtonian fluid flowing in a cylindrical tube, which mimics 
a natural pore, whose walls are oscillating longitudinally and the fluid is subject to an 
oscillatory pressure gradient. We give analytical solutions of the problem in the frequency 
domain. 

The governing equation of the problem is the linearized momentum equation of an in- 
compressive fluid 

/)- = -Vp + /iVV (1) 

Here, v, p, p denote velocity, pressure and mass density of the fluid, while \x is the viscosity 
coefficient. 

Now, let u be a velocity of the wall of the tube which oscillates in time as e~"* w *. The 
flow of fluid in a cylindrical tube with longitudinally oscillating walls can be described by a 
singe component of the velocity, namely, its ^-component v z (z axis is along the centerline 
of the tube). We use cylindrical coordinate system (r, 0, z) in treatment of the problem. We 
introduce the relative velocity U\ as U\ — v z — u. Thus, assuming that all physical quantities 
vary in time as e~ tujt , we arrive at the following master equation for U\ 

V ^ + "Lu, = -1. (2) 
v v 

Here, we have introduced the following notations: 

pX = -(Vp + p— ), 



which is a sum of the applied pressure gradient and force exerted on the fluid from the 
oscillating wall of the tube and, u, which is v — fi/p. 
The solution of Eq.(2) can be found to be ]?) 

U 1 (r) = ~ + CJ (Pr), 
too 



where Jo is the Bessel function and (3 = Jiuj/v. 

Assuming that the slip velocity is U\(a) = U s at the wall of the tube, where a is its radius, 
we obtain 



U x (r) = -- 
iuj 



1-(1 + C7)W 
[l + Us) J ((3a) 



(3) 



Here. 



U„ = U^ = U, 



X ~°a 2 X 
Defining the cross-section averaged velocity as 

2 r a 



((3a) 2 . 



2 f a 

Ui = — Ui(r)rdr, 
cr jo 



we obtain 



Xa 2 



2(1 + U s )J 1 ((3a) 
(J3a)J (J3a) 



v ((3a) 2 

Following work of Biot |7|] we calculate the stress at the wall r, 

'dU^r^ n(3X fi , T - \ Ji(/3a) 



T = —jl 



(l + Us) 



(4) 



(5) 



, 9r J r=a iuj v J ((3a)' 
The total friction force is 2-nar. Following Biot we calculate the ratio of total friction 
force to the average velocity, i.e. 

2wafx o _ fo N i r - ^ J\((3a) 



x 



-2nfi((3a) (l + U s ) 



2(1 + U S )J 1 ((3a 



J ((3a) 



-i 



(6) 



((3a)J ((3a) 

Simple analysis reveals that (assuming U s — - > as uj — > 0, see discussion below) 

2-aar 
hm — = — = 87m, 

which corresponds to the limiting case of Poiseuille flow. Following Biot |7j , we also introduce 



a function F(k) with k being frequency parameter, k = aycu/V, in the following manner 

2naT 



U-i 



M^F(k) 



thus, 



F(k) 



1 



Vi (l + U a ) 




2J 1 (K V / i) 



n -1 



X 1 




K\/lJo(K\/l) 



(7) 



Note, that F(k) measures the deviation from Poiseuille flow friction as a function of fre- 
quency parameter k. The Biot's expression for F(k) in the no boundary slip regime can be 
easily recovered from Eq.(7) by putting U s — > for all k's. 

So far, we did not specify U s , however there are certain physical constraints it should 
satisfy: 

(A) Authors of Ref. |1|] demonstrated that the slip length (which, in fact, is proportional 
to the slip velocity) is a function of the approach rate and they showed why several previ- 
ous careful measurement of confined liquids have not observed evidence for boundary slip. 
Under the low approach rates employed in previous measurements slip length is zero and 
no-slip boundary condition is applicable. Experiments reported in Ref. [||] were performed 
when half-sphere approached a plane with liquid placed in-between at different approach 
rates. However, we should clearly realize what term "approach rate" means in the context 
of Biot's theory: Biot investigated fluid flow in the cylindrical tube whose walls are har- 
monically oscillating in the longitudinal direction as x(t) = Ae~ lwt , therefore if similar to 
Ref. m experiment would be done for the oscillation tube, the " approach rate" would be the 
amplitude of x(t), i.e. —iu>A. Thus, when uj — > 0, U s should also tend to zero. 

(B) At high frequencies the viscous effects are negligible compared to the inertial effects. 
Thus, fluid at high frequencies behaves as an ideal (non-viscous) pore fluid, which allows 
for an arbitrary slip. Therefore, when uj — > oo, U s should tend to zero. When we examine 
Fig. 6 from Ref. pj and Fig. 3 (bottom panel) from Ref. closely we should understand that 
the authors measured ascending part where slip length increases with the increase of the 
approach (flow) rate. Since, the viscous effects should be negligible compared to the inertial 
effects for large approach rates, there must be also a descending part - clearly slip length (or 
the same as slip velocity) cannot grow infinitely as approach rate increases. Viscous effects 
should give in to the inertial effects and fluid should behave an ideal fluid allowing for an 
arbitrary slip. 

Based upon the above two physical arguments we conclude that U s should be a smooth 
function which tends to zero when uj (or the same as n) tends to both zero and oo. Therefore, 
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we postulate following phenomenological expression for U s 

Re[U s {K)] = Im[U s {K)}=i^^- 



2^5 



This function has a maximum at = y2A/3. By fixing B at 5 5 v4 3 /(4 x 3 3 ) we force it 
to achieve its maximum equal to unity at k = k» (when £ = 1). We introduced £ (usually 
< £ < 1) as a sort of "weight" of the boundary slip effect into the solution. Putting, £ = 
this way would easily allow us to recover no-slip case, while the increase of £ we would be 
able to trace effect of non-zero slip onto the all physical quantities. We plot this function in 
Fig. 1 for A = 25 and f = 1. 

It is worthwhile to mention that, at first glance, it seems that non-zero boundary slip, 
which appears at "high approach rates", should be attributed to the non-linear effects. 
However, thorough interpretation of the experimental results of Ref. in the context of 
the oscillatory tube (see points (A) and (B) above) allows us to conclude that the non-zero 
boundary slip can be also incorporated into the linear Biot's theory. 

Next, we plot both F t (k) = Re[F(n)] and F^k) = —Im[F(K)] in Fig. 2 for the three 
cases: when there is no boundary slip (£ = 0), and £ = 0.05, 0.1. We gather from the plot 
that the £ = is identical to Fig. 4 from Ref.]7) as it should be. However, we also notice a 
noticeable difference from the classic case when £ is non-zero in the intermediate frequencies 
domain. Of course, according to our definition of the phenomenological form of U s (k), even 
for non-zero £ (non-zero boundary slip), when k — > and k — > oo, asymptotically F(k) 
behaves as classic Biot's solution 0, i.e. 



limF r (K) = 1, (8) 



and 



limF i (K) = 0, (9) 

^F( K ) = -Vi=- i (-^-y (io) 

Since in our phenomenological model we allow for the deviations in the intermediate fre- 
quency domain, it is easy to foresee that these will have an impact on all of the predictions 
of the Biot's theory precisely in that frequency range. Namely, all observable quantities 
predicted by the Biot's theory, such as phase velocities and attenuation coefficients of the 
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both rotational and dilatational waves will be affected by the introduction of boundary slip 
into the model in the intermediate frequency range. 

Biot showed that the general equations which govern propagation of rotational and 
dilatational high-frequency waves in a fluid-saturated porous medium are the same as in the 
low-frequency range provided the viscosity is replaced by its effective value as a function of 
frequency. In practice, it means replacing the resistance coefficient b by bF(n). 

The equations describing dynamics of the rotational waves are 

— (p u u + Pl2 ti) + bF(K)-(uj -fi) = NV 2 u, (11) 

-^(p 12 ^ + P22^)-bF(K) — (uj-n) = 0, (12) 

where, pm P12 an d P22 are mass density parameters for the solid and fluid and their inertia 
coupling; uj = curl u and Q = curl U describe rotations of solid and fluid with u and U being 
their displacement vectors, while the rigidity of the solid is represented by the modulus N. 
Substitution of a plane rotational wave of the form 

uj = de i{lx+xt \ Q = C 2 e i{lx+Xt \ (13) 

into Eqs.(ll) and (12) allows us to obtain a characteristic equation 

Nl 2 



pa 2 



Er-iEi, (14) 



where I is wavenumber, \ — is wave cyclic frequency, p = pn + %Pi2 + P22 is the mass 
density of the bulk material and a is a pore radius. 

The real and imaginary parts of Eq.(14) can be written as 

„ (711722 - 7i 2 2 ) (722 + £2) + 72262 + 4 + 4 f , K x 

(722 + e 2 ) 2 + e( 



and 

£1(712 + 722) 2 
(722 + e 2 ) 2 + ef 

where 7^ = pij/p, e x = (712 + 722) (/c//) F t (k) = (712 + 722) (/ c //) F r (5yjf/f c ), e 2 = (712 + 



p tU7i2 -r 122) 



722) (/ c / 7) -^( K ) = (712 + 722) (/c//) Fi(5y/f/f c ). The function was written here more 
conveniently as a function of frequency / , i.e. F(k) = F(SJ f / f c ) 0, where 5 is a factor 
dependent on pore geometry. For the hollow cylinder-like pores, S = a/8 [0] and we use 
this value throughout the paper. f c is the critical frequency above which the Poiseuille flow 



breaks down, and it equals b/ (2iip 2 ). Here p 2 denotes the product of porosity and fluid mass 
density. 

In order to obtain phase velocity and attenuation coefficient of the rotational waves, we 
put I = Re[l] + ilm[l]. Thus, the phase velocity is then v r = x/|-R e K]|- Introducing a 



reference velocity as V r = JN/p, we obtain the dimensionless phase velocity as 

V ~ r ; 1 lV2" 

! E? + El + E r 



(17) 



To obtain the attenuation coefficient of the rotational waves, we introduce a reference 
length, L r , defined as L r = V r /(2-Kf c ). The length x a represents the distance over which the 
rotational wave amplitude is attenuated by a factor of 1/e. Therefore we can construct the 
dimensionless attenuation coefficient as L r /x a , 

lV2 



Ef + El - E r 



X a fc \/2 



(18) 



The equations describing dynamics of the dilatational waves are 



V 2 (Pe + Qe) = — ( Pll e + Pl2 e) + bF{K)-{e - e), (19) 

V 2 (Qe + Re) = ^(p 12 e + p 22 e) - bF{K)-{e - e), (20) 

where, P, Q and R are the elastic coefficients, e = div u and e = div U are the divergence of 
solid and fluid displacements. Again, substitution of a plane dilatational wave of the form 

e = C ie i{lx+Xt \ e = C 2 e i{lx+Xt \ (21) 

into Eqs.(19) and (20) allows us to obtain a characteristic equation 

(z - Zl ) (z - z 2 ) + iM(z - 1) = 0, (22) 

where z = l 2 V^/x 2 , V 2 = (P + R + 2Q) / p represents the velocity of a dilatational wave 
when the relative motion between fluid and solid is absent, = V 2 /V 2 2 with V\,i being the 
velocities of the purely elastic waves with subscripts 1,2 referring to the two roots of Eq.(22), 
and finally M = (e 1 + ie 2 ) / \a u a 22 - aj 2 ) with a n = P/(P + R + 2Q), a 22 = R/(P + R + 2Q) 
and (7i2 = Q/(P + R + 2Q). 
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Eq.(22) has two complex roots zj and zu. Phase velocities of the two kinds of dilatational 
waves can be defined as 



vj 1 vn 1 



V c Re[^zJY V c Re[^/zHY 
while the corresponding attenuation coefficients can be also introduced as 



(23) 



^ = Im\^\L ^L = /m[v /i77]Z. (24) 

Xl Jc Xll Jc 

III. NUMERICAL RESULTS 

In order to investigate the novelties brought about into classical Biot's theory of prop- 
agation of elastic waves in porous medium by the inclusion of boundary slip, we have 
studied the full parameter space of the problem. 

In all forthcoming results, we calculate phase velocities and attenuation coefficients for 
the case 1 from Table I taken from Ref. f7|, which is o~n = 0.610, 022 = 0.305, cr 12 = 0.043, 
7ll = 0.500, 722 = 0.500, 712 = 0, z x = 0.812, and z 2 = 1.674. 

We calculated normalized phase velocity of the plane rotational waves, v r /V r , and the 
attenuation coefficient L r /x a using our more general expression for F(k) (which takes into 
account non-zero boundary slip) given by Eq.(7). 

In Fig. 3 we plot phase velocity v r /V r as a function of frequency for the three cases: the 
solid curve corresponds to £ = (no boundary slip), while long-dashed and short-dashed 
curves correspond to £ = 1.0 and £ = 1.5 respectively (we have used these large values in 
order to emphasize the effect of the boundary slip). Note that the £ = case perfectly 
reproduces the curve 1 in Fig. 5 from Ref. |7[]. For £ = 1.0 and 1.5 we notice a deviation 
from the classic behavior in the form of an decrease of phase velocity in the domain of 
intermediate frequencies. 

Fig. 4 shows the attenuation coefficient L r /x a of the rotational wave as a function of 
frequency for the three values of £: the solid curve corresponds to £ = (no boundary 
slip), while long-dashed and short-dashed curves represent the cases £ = 1.0 and £ = 1.5 
respectively. Note that £ = case coincides with curve 1 in Fig. 6 from Ref. |7j. We observe 
that, in the intermediate frequency range, with increase of boundary slip our model yields 
substantially higher values of the attenuation coefficient than in the classic no-slip case, 
which is in qualitative agreement with the experimental data ||, H|. 
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We also calculated normalized phase velocities of the plane dilatational waves, Vj/V c and 
vii/Vc, and the attenuation coefficients L c /xi and L c /xn using our more general expression 
for F(k) given by Eq. (7). 

In Figs. 5 and 6 are similar to Fig. 3 except that now we plot phase velocities vj/V c and 
vujV c as a function of frequency for £ = 0, 0.05, 0.1. Note that the solid curves on the both 
graphs reproduce curves 1 in Figs. 11 and 12 from Ref. ]7[|. We gather from the graph that 
introduction of the non-zero boundary slip leads to decrease of vi/V c , while Vn/V c increases 
in the in the domain of intermediate frequencies as boundary slip, £, increases. 

In Figs. 7 and 8 we plot the attenuation coefficients L c j xj and L c j xn in a similar manner 
as in Fig. 4, but now for £ = 0, 0.05, 0.1. Again we observe that the £ = case reproduces 
curves 1 in Figs. 13 and 14 from Ref.j7j. We gather from Figs. 7 and 8 that introduction 
of the non-zero boundary slip yields higher than in the no boundary slip case values of the 
attenuation coefficients for both types of the dilatational waves in the domain of intermediate 
frequencies, which, again, is in qualitative agreement with the experimental results f|, @. 

It is worthwhile to note that for no-zero £ the asymptotic behavior of the elastic waves 
in the f / f c — * oo as well as f / f c — > limit is identical to the classic behavior established by 
Biot J7]. This is a consequence of the assumptions of our phenomenological model, which is 
based on robust physical arguments. 

IV. DISCUSSION 

In this paper we have studied the effect of introduction of the boundary slip in the theory 
of propagation of elastic waves in a fluid-saturated porous solid originally formulated by 
Biot [0]. Biot's theory does not account for boundary slip effect, however, the boundary slip 
becomes important when the length scale over which the fluid velocity changes approaches 
the slip length, i.e. when the fluid is highly confined, for instance, fluid flow through porous 
rock or blood vessel capillaries. In the light of recent convincing experimental evidence of a 
boundary slip in a Newtonian liquid 0, it is necessary to take into account this effect into 
the Biot's theory where appropriate. We have studied the effect of introduction of boundary 
slip upon the function F(k) that measures the deviation from Poiseuille flow friction as a 
function of frequency parameter k. This function crucially enters Biot's equations which 
describe dynamics of fluid-saturated porous solid. Therefore, a revision of Biot's theory 
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was needed in order to incorporate boundary slip effect into the all measurable predictions 
of this theory such as phase velocities and attenuation coefficients of the both rotational 
and dilatational waves. We have performed such analysis, and in summary, we found that 
the introduction of the non-zero boundary slip into the Biot's theory of propagation of the 
elastic waves in a fluid-saturated porous solid results in 

• the decrease, as compared to the no-slip limiting case, of the phase velocities of both 
rotational waves (v r /V r ) and dilatational wave of the first kind (vi/V c ) in the domain 
of intermediate frequencies. On contrary, the phase velocity of the dilatational wave 
of the second kind (vu/V c ) experiences an increase as compared to the no-slip limiting 
case in the domain of intermediate frequencies. 

• in the domain of intermediate frequencies the attenuation coefficients of both the ro- 
tational (L r /x a ) and dilatational waves (L c /xj and L c /xn) are increased as compared 
to the no-slip limiting the boundary slip increases, which is in qualitative 
agreement with the experimental data |§, @. 

• behavior of all physical quantities which describe the elastic waves in the asymptotic 
limits of both small and large frequencies is not affected by the introduction of the 
non-zero boundary slip. The deviation occurs only in the domain of intermediate 
frequencies, as prescribed by our phenomenological model. 

The investigation of properties of elastic waves is important for a number of applications. 
The knowledge of phase velocities and attenuation coefficients of elastic waves is necessary, 
for example, to guide the oil-field exploration applications, acoustic stimulation of oil pro- 
ducing fields (in order to increase the amount of recovered residual oil), and the acoustic 
clean up of contaminated aquifers || |], ||, Therefore, our results would be useful for 
various applications in oil production as well as in ecology. 

From the recent experimental results of Ref.|2| we gathered that there are physical sit- 
uations were the no-slip boundary condition becomes invalid. We have formulated a phe- 
nomenological model of elastic waves in the fluid-saturated porous medium based on Biot's 
linear theory and certain physically justified assumptions on the variation of boundary slip 
velocity with frequency, U s (k). Since, there are no experimental measurements of U s (k), for 
a cylindrical tube, on which "Biot-like" theory relies, there is a certain freedom of choice, 
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which could be used to obtain a better fit of experimental data with the theory in cases 
where classic Biot's theory fails to do so. If fact, our model predicts higher than the Biot's 
theory values of attenuation coefficients of the both rotational and dilatational waves in the 
intermediate frequency domain, which is in qualitative agreement with the experimental 
data 0. Therefore, the introduction of the boundary slip yields three-fold benefits: 

• Better agreement of theory with an experimental data, since, the parametric space of 
the model is larger (includes effects of boundary slip). 

• Possibility to identify types of porous medium and physical situations where boundary 
slip is important. 

• Constrain model parameters that are related to the boundary slip. 

We would like to close this paper with the following two remarks: First, the reported 
slip length varies from tens of nano- meters |L| to a few microns JJ, depending on a type of 
a liquid and differences in the experimental set up. The latter experimental work clearly 
demonstrates the applicability of our model to the elastic wave phenomena in a usual porous 
media found in Nature. Second, in this work we postulated a phenomenological expression 
for the boundary slip velocity based on some physically sound arguments (in fact, similar, 
successful, approach has been used by the authors of Ref.|To| in the context of solar physics). 



However, perhaps it may be possible to derive an expression based on more general "first 
principles" , which should be a subject of a future work. 
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Figure Captions 

Fig. 1: Behavior of Re[U s (n)] = Im[U s (K)] as function of frequency parameter, k, for 
f = 1, A = 25. 

Fig. 2: Behavior of F r (n) (thick curves) F^k) (thin curves) as function of frequency 
parameter, k, according to Eq.(7). Solid curves correspond to the case when £ = (no 
boundary slip), while long-dashed and short-dashed curves correspond to £ = 0.05 and 
£ = 0.1 respectively. 

Fig. 3: Behavior of dimensionless, normalized phase velocity of the rotational wave, v r /V r , 
as a function of frequency. Solid curve corresponds to the case when £ = (no boundary slip), 
while long-dashed and short-dashed curves correspond to £ = 1.0 and £ = 1.5 respectively. 

Fig. 4: Behavior of dimensionless, normalized attenuation coefficient of the rotational 
wave, L r /x a , as a function of frequency. Solid curve corresponds to the case when £ = 
(no boundary slip), while long-dashed and short-dashed curves correspond to £ = 1.0 and 
£ = 1.5 respectively. 

Fig. 5: Behavior of dimensionless, normalized phase velocity of the dilatational wave, 
vi/V c , as a function of frequency. Solid curve corresponds to the case when £ = (no 
boundary slip), while long-dashed and short-dashed curves correspond to £ = 0.05 and 
£ = 0.1 respectively. 

Fig. 6: The same as in Fig. 5 except for the curves are for v n /V c . 

Fig. 7: Behavior of dimensionless, normalized attenuation coefficient of the dilatational 
wave, L c /xi, as a function of frequency. Solid curve corresponds to the case when £ = 
(no boundary slip), while long-dashed and short-dashed curves correspond to £ = 0.05 and 
£ = 0.1 respectively. 

Fig. 8: The same as in Fig. 7 except for the curves are for L c jxu. 
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